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p^, ABSTRACT 
O 

' A general physical mechanism of the formation of line-driven winds at the 

■ vicinity of strong gravitational field sources is investigated in the frame of General 
Relativity. We argue that gravitational redshifting should be taken into account 

^ ' to model such outflows. The generalization of the Sobolev approximation in the 

■ frame of General Relativity is presented. We consider all processes in the metric 
of a nonrotating (Schwarzschild) black hole. The radiation force that is due 
to absorbtion of the radiation flux in lines is derived. It is demonstrated that 
if gravitational redshifting is taken into account, the radiation force becomes 
a function of the local velocity gradient (as in the standard line-driven wind 
theory) and the gradient of goo- We derive a general relativistic equation of 
motion describing such flow. A solution of the equation of motion is obtained and 
confronted with that obtained from the Castor, Abbott & Klein (CAK) theory. 
It is shown that the proposed mechanism could have an important contribution 
to the formation of line-driven outflows from compact objects. 



Subject headings: radiation mechanisms: general - stars: mass loss - stars: winds, 
outflows - galaxies: active 
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Introduction 

It has been demonstrated both theoretically and observationally that accretion disks 
around compact objects can be powerful sources of fast plasma outflows. Among the most 
important processes known to work are magnetic and radiation driving. In fact radiation- 
driven winds can exist in most of systems where accretion disk can produce enough UV radi- 
ation (the standard line-driven wind theory gives an approximate value of Luv > lO^^i^edd). 
However, this conclusion should be treated with care because of the physical conditions in 
a disk wind which are very different from that of an 0-typc star wind. While realistic ac- 
cretion disk winds are most likely driven by the combination of the radiation and magnetic 
forces here we focus on the scenario when momentum is extracted most efficiently due to 
absorption of the radiation flux in lines of abundant elements. 

In the paper of Dorodnitsyn (2003) (hereafter Dl) it was proposed a mechanism when 
line-driven acceleration occurs in the vicinity of compact object so that the the gravitational 
redshifting can play an important role. The generalization of these studies in the frame of 
General Relativity (GR) is the problem that we address in this paper. A mechanism that 
we study is quite general and can be considered to work in any case when there is enough 
radiation to accelerate plasma and radiation driving occurs in strong gravitational field. 
Particularly we discuss winds in active galactic nuclei as they manifests most important 
properties of accretion disk -|- wind systems keeping in mind however that our treatment 
allows to consider their low mass counterparts. 

It is widely accepted that a supermassive black hole (BH) hes in the cores of most of 
active galactic nuclei (AGN). The accretion activity around such a black hole results in a 
production of a powerful continuum radiation - a defining characteristic feature of the quasar 
phenomenon. The dynamical role of this radiation is so high that it is probably responsible 
for the formation of fast winds which are observed in AGN. The radiation pressure on lines 
plays the crucial role in acceleration of such outflows. The most prominent feature seen in 
about 10% of quasar spectra are the broad absorbtion line systems BALs - the blue-shifted 
UV resonance lines from highly ionized species (N V, CIV, Si IV). These come from ions of 
differing excitation with bulk velocities of up to 0.2 c. A successful model must also explain 
a simultaneous existence of NALs - narrow absorbtion line systems (N V, CIV), seen in UV 
and X-rays from about a half of Seyfert galaxies and associated with outfiows of 1000 kms~^, 
and BEL - broad emission lines present in all AGNs indicating fiows as fast as 5 • lO^kms"^. 
These well established features together with total luminosity of up to L ~ 10^^ ergs~^ gives 
us the crucial evidence of the dynamical importance of the line-driven mechanism in AGNs. 
A quasi-lD model of the quasar wind was developed in Murray et al. (1995) (however it 
is not clear how justified is the assumption that equations in radial and polar directions 
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could be solved separately). The 2D calculations of the accretion disk powered winds were 
made in Proga, Stone & Drew (1998) and Proga, Stone & Kallman (2000), while the winds 
from massive X-ray binaries (together with ionizing effects of the radiation from the central 
source) have been considered in Stevens & Kallman (1990). 

In the pioneering paper of Sobolev (1960), it was recognized that the problem of the 
radiation transfer in lines in a continuously accelerating medium is simplified drastically 
in comparison with that of a static case. In the paper by Lucy & Solomon (1970) it was 
pointed out the importance of the line opacity for the formation of winds from hot stars. 
Most of our understanding of the line-driven mechanism is based on the prominent paper of 
Castor, Abbott and Klein (1975), (hereafter CAK) where a theory of the 0-type star wind 
was developed. These studies explained how a star that radiates only a tiny fraction of its 
Eddington limit, can have a very strong wind. CAK was able to demonstrate that the radi- 
ation force from en ensemble of optically thin and optically thick lines can be parameterized 
in terms of the local velocity gradient. This elegant theory was further developed in papers 
of many authors. All this work resulted in what is usually called a "standard line-driven 
wind theory" (hereafter SLDW). 

It is rather problematic, however, to apply directly the CAK theory to accretion disk 
winds because of the geometrical difference and because of the different properties of the 
spectrum emitted by the central source of the continuum radiation. For example, a wind in 
AGN is likely to be launched from accretion disk-like structure and thus is intrinsically two 
dimensional with the geometry that is close to axial symmetry. The second crucial difference 
is that in active galactic nuclei a wind is exposed to a hard UV and X-ray continuum radiation 
that stripes electrons from abundant elements much more effectively than the quasi-black- 
body radiation of the hot luminous star. In case of AGN the radiation flux produces highly 
ionized species over much of the wind. The considerable lack of the atomic data for highly 
ionized ionic species and of intrinsically 2D radiation-hydro and transfer calculations makes 
a task of the reahstic modeUing of AGNs winds very problematic. 

In the standard line-driven wind theory a given parcel of gas sees the matter that is up- 
stream redshifted because of the difference in velocities (assuming that a wind is accelerating 
gradually). This helps a line to shift from the shadow produced by the underlying matter 
and to expose itself to the unattenuated continuum. It was shown in Dl, that together 
with Sobolev effect the gravitational redshifting of the photon's frequency should be taken 
into account when calculating the radiation force. In case of strong gravitational field the 
gradient of the gravitational potential works in the same fashion as the velocity gradient 
does when only Sobolev effect is taken into account, so that the radiation force becomes 
gi ~ {dv/dr + ^d(j)/dr). As it was shown in Dl now the gravitational field works in ex- 
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posing the wind to unattenuated radiation of the central source. Thus we call such a flow 
"Gravitationally Exposed Flow" (GEF). 

Conditions present in the inner parts of the realistic accretion disk wind are far from 
being clear. However it is known that most of the radiation flux is produced in the innermost 
parts of the accretion disk. We may expect that some part of the wind which is located 
beyond few tens of gravitational radii may be moving quasi-radially. To make our treatment 
as general as possible, we consider spherically - symmetrical radiationally accelerated wind. 
Lines and electron scattering are assumed to be the only sources of opacity, no ionization 
balance is calculated. 

In Dl gravitational field was considered by means of the gravitational potential. Thus 
all equations were in fact derived in the fiat space-time, and only the effect of the gravi- 
tational redshifting was taken into account when calculating the radiation pressure term. 
The resultant solution was then compared with CAK wind solution. This approach is not 
self-consistent. In GR when calculating the radiation force, the effect that is due to Doppler 
shifting should be taken into account simultaneously with gravitational shifting (no bending 
of photon trajectories is considered since the force is calculated in radial streaming limit). 
Obviously, the CAK - type solution can exist only in the flat spacetime. It is important to 
note that the self- consistent modelling of GEF is possible only via general relativistic treat- 
ment. Thus, the main goal of this paper is to compare the general relativistic GEF solution 
with the SLDW solution, obtained in the Newtonian gravity. 

Here we solve GR equations of motion for radiatively accelerated wind and calculate 
the radiation pressure in the radial streaming limit in the Sobolev approximation. Making 
use of the Sobolev approximation allows us not to treat the General Relativistic radiative 
transfer formalism. There exist an extensive literature, where the radiative transfer problem 
in GR is considered for the purpose of hydrodynamical calculations of spherically symmetric 
accretion (e.g. TuroUa & Nobili 1988; Thorne, Flammang Sz Zytkow 1981; Nobili, TuroUa, 
& Zampieri 1991). In these papers the radiative moment formahsm of Thorne (1981) had 
been extensively used. However, for our purposes it is not needed to use this sophisticated 
formalism. In the Sobolev approximation the flow is treated in fact as non-relativistic and 
the only source of the opacity is the line and electron scattering. In such an approach it is 
possible to derive the radiation force without an explicit solution of the radiative transfer 
equation. Thus we use only escape probability arguments - exactly as the radiation force 
is derived in the standard line-driven wind theory (e.g. Mihalas 1978; Lamers & Cassinelli 
1999) 

The plan of this paper is as follows. In Section 1 the General relativistic equations, 
describing the matter that is interacting with the radiation field are derived. Then, in 



- 5 - 



Section 2 the optical depth and the radiation force in Sobolev approximation are obtained. 
We obtain the equation of motion for the hne-driven wind in GR and then numerically 
calculate a set of its solutions in Section 3. The results are summarized and the perspectives 
are discussed in the conclusions. 



1. The radiation-driven wind 

It was suggested in Dl that taking into account the gravitational redshifting in mod- 
elling the line-driven winds can substantially increase the efficiency of the mechanism. 

The wind is assumed to blow in the background metric described by the familiar 
Schwarzschild line element 

ds'^-hdi' + ^ + r'^dn'', (1) 

where h= (l — ^) , = 2M, where M is a black hole mass. Geometrized units (c — G — l) 
are used throughout this Section. Let O be an observer that is at rest at r, that measures 
physical quantities like -0- velocity, dl = dr/Vh - displacement etc. The observer 6 has 
the following tetrad of orthonormal basis vectors: e£ = /i~^^^et, if = h^^'^ er, = ee, 
= (r sin 6')"^ 60. 

The stress-energy tensor for the ideal gas reads: 

rpaP ^ (p + p-^u°'u'^ + PC/"'^, (2) 

where for the total mass-energy density we have 

p^Po{l + Ei), (3) 

where po = nT-bUb is the barionic rest-mass density, is the baryon number density, and Ei 
is the internal energy of the ideal gas per unit mass. 

The continuity equation reads: 

{nbu"),a = 0. (4) 

The process of the interaction of the matter with the radiation is described by the four-force 
density which is given by 

G" = y dv j dVL{xJv-r]u)n'', (5) 
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where Xv is the opacity and r]v is the emissivity. 1^ is the specific intensity of photons of the 
frequency propagating in the direction n", the space component gives the net rate of 
the radiative momentum input, while cG° equals the rate of the radiative energy input. 

The equations of hydrodynamics for the matter interacting with the radiation filed read 
(see, e.g., Mihalas & Mihalas, 1984): 

Tf - G"* = 0. (6) 

Applying the projection tensor 

to the equation (6) will result in the Euler equation of motion: 

Pa,{T°'f-G'^]^Q. (8) 
Calculating equation (8) in the orthonormal frame of O we obtain: 

P + pf^^dv M\ dP 1 ^5, ^ 

' ^t'^^I + 117 + ^ - - ^G') = 0, (9) 



h \ dr J dr 

where 7 = (1 — w^)"^/^ and it was taken into account that the radiation force measured by 
O is given by: = G^/\/h, G° = VhG^. These physical components of the radiation 
four-force density can be further transformed to the frame that instantly coincides with the 
moving gas by making use of the corresponding Lorentz transformation: j{G^ — v G^) — Gl- 
The relative importance of the frame-dependent term G^ in (9) is 0{v/c) relative to G^ 
(Mihalas & Mihalas 1984). The radiation force G^ = G^ ^ + Gq is the sum of a radiation 
force due to absorbtion in continuum and a force which results from the line transition, as 
measured by the physical co-moving observer. The latter should be calculated taking into 
account both the shifting of frequency due to Dopplcr effect and the gravitational redshifting. 
Note that equation (9) corresponds to equation [11] of Nobili, TuroUa, & Zampieri (1991), 
if to imply that G^ consists only of the part of the force that is due to continuum radiation 
fiux. Prom equation (4) we obtain continuity equation in the form: 

Povr^-fVh = — , (10) 

where (7P0) - is the barionic rest-mass density as measured in the laboratory frame. 

The key ingredient of the CAK theory is that the optical depth which is due to the fine 
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absorbtion can be expressed as a function of the velocity gradient in the flow. The same is 
true in the case when a photon, emitted somewhere deeply in the potential well, will become 
resonant with the line absorption both due to velocity difference and due to GR effect of 
redshifting. 



2. Optical depth and radiation force in Sobolev approximation 

A photon emitted at a given radius will suffer a continuous both gravitational and 
Doppler redshifting and may become resonant with a line transition at some point down- 
stream. Thus a ray of a frequency f^, emitted by the matter that for simplicity is assumed 
to be at rest at radius r^^, at a given point r has a frequency Uiab, as measured by the observer 
O and that is obtained from relation (see, e.g. Landau & Lifshitz 1960): 



^dV l^d — ^lab 



VW)-^"^. (11) 



where y°° is the frequency of the ray at infinity. We restrict ourself to the radially streaming 
photons only and assume that they arc emitted from a point source. In such a case the 
Sobolev optical depth may be calculated without solving the radiation transfer equation. 
Optical depth between and a given point r can be written (Novikov & Thorne 1973): 

r - r dr 

= / Xl,labdl= / Xl,lab , _ 2MN1/2 ' (-'-^) 

where dl is the proper-length element, and Xi,iab (cm~^) is the absorbtion coefficient as mea- 
sured by O. However, it is more appropriate to measure absorbtion (as well as the emission) 
in the local frame co-moving with the fluid. In such a case, opacity is transformed according 
to the reation: Xi,iab = Xi,com ^/^lab, where in the co-moving frame the frequency of the ray is 
u — 7i/;ab(l — /?), /3 = v/c, and Xi,com is the absorbtion coefficient measured by the co-moving 
observer. In the co-moving frame, the line-center opacity is determined by the following 
relation: 

iV./9.-iV./.. _ (13) 

mc Aud 

where Aud = ^oVth/c is the Doppler width, and uq is the line frequency, / is the oscillator 
strength of the transition, g is the statistical weight of the state, A^"^, and gu, gi are 
respective populatios and statistical weights of the corresponding levels of the line transition. 
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Then, one have to compare 9 with the frequency of the hne uq, to see whether it is in the 
range of the hne profile (p{iy — i/q). That is in fact a standard procedure that is followed in 
order to calculate the Sobolev optical depth, but with an additional step - the gravitational 
redshifting. Thus, both Doppler and gravitational redshifting of the photons frequency are 
taken into account. 

Introducing a new frequency variable: 



y = i>-iyo^-f{l- P)-j=-iyo, (14) 
We change the integration variable in (12) from r to y: 



ip{iy - i^oJxA^D ^ = T-i ^dy, (15 

4 ''lab Vh Jyira) ''lab ^oo ( ^ _ 

where 7] = 7(1 — f3). The relation (14) was used in order to calculate dx/dr. The last term 

in the denominator of (15) is due to the "gradient of the gravitational field": w = —\fh. 

GM 

Note that w • = ^ - equals the acceleration of the free-faUing particle that was initially 

at rest in the Schwarzschild metric. In case of w = we obtain the result of Hutsemekers 

y^Vthi^ — S) 

& Surdej (1990): ti — — — — — . Assuming that the hne profile is a 8- function, or, 

^dvjdr 

equivalently, that the region of interaction is infinitely narrow, we find the optical depth 
in Sobolev approximation: = — , where it was taken into account that 



Vh'-f— + W7 
dr 

f]' = —7^(1 — 13)13' . In our treatment we retain only terms 0{v/c) (in the equation of motion) 
and thus resultant Sobolev optical depth can be written in the form: 



r-dv 

V h— + cw 
dr 

where (cm^-g~^) is the mass absorbtion coefficient: k'i = x?/PO) measured in the rest-frame 
of the fiuid, and 7 = 1. 

In the weak field limit the optical depth (16) will transform to equation (11) of Dl. If 
there is no gravitational redshifting taken into account then the Sobolev optical depth is ob- 
tained: Tsob = I = I'sob, and the Sobolev length scale Isob = Vth/ (dv/dr ) determines 
\dv/dr\ 
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a typical length on which a line is shifted on about its thermal width. In general case from 
(16) we conclude that 



where Iqef 



n^x'ihEF, (17) 

Vth 



r-dv 

Vh— + cw 

dr 



The radiation force form a single line exerted by the material as measured in its rest 
frame reads: 

5. ^ f.°(-)x?/p.A^.>l-e-_ 
c n 

where F°(i/;)(erg • cm~^ • Hz~^ • s~^) is the radiation flux at the line frequency in the rest- 
frame of the fluid. Note that ti in (18) should be calculated with taking into account the 
redshifting and is given by (16). The e""^' term reflects the fact that the incident flux at r is 
reduced in comparison with the initial flux Fj,. 1 — e""^' gives the "penetration probability" 
for a ray to reach a given point. 

In our treatment we neglect special relativistic terms all equations. We may expect 
that final results will be at least qualitatively correct for the flows as fast as ~ 0.2 0.3c. 
Prom (18), we conclude that if a line is optically thin, the radiation force does not depend 

upon the redshifting law (11): gmn — — —Fi,. On the other hand, when t; >> 1 the 

optically thick line produces force that can be roughly described by the following: g^thick = 

— T~^gihin, and thus it is independent of the line strength and gi ~ dv/dr. 

cpo r 

The difference between po,Za6 and po is 0{v'^/c^). Note, that in our case g^^^^ — x°F°/c = 
5'thin + 0{v'^/c'^), 5(°hick = Fl^/c = giuck + 0{v'^/c^), and a factor (1 - (3) was omitted when 
calculating (16). 

According to the well accepted notation let introduce the optical depth parameter: 

CTe PoVth 



r-dv 

V h— + cw 

dr 



(19) 



where is the electron scattering opacity per unit mass, and t is connected to r/ via relation: 
Ti = C,t, where = K^ae- 

The role of the parameter ^ is very important because it allows to separate the line 
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optical depth into two parts: the first (^) that depends on statistical equihbrium, and the 
second (t) which depends only on the redshifting law (14). It will alow us to use the standard 
parameterization law for the force multiplier when calculating the radiation force. Summing 
(18) over the ensemble of of optically thin and optically thick lines we obtain total radiation 
acceleration: 

9i-J29i-^Mit)^ (20) 

I 

where F is the total flux and the "force multiplier" M(t) equals: 

Ti<l Tl>l 

CAK found that M(t) can be fitted by the power law: 

M{t) = kr''. (22) 



The equation of motion 



Combining equations (19), (20), (22) and the equation of continuity (10), the equation 
of motion describing a stationary, spherically-symmetric, wind can be cast in the form: 



GM 



dr ' r'^ J po df^ -/hAnr'^c 
y/h4nr^c \avthM J 



Vhvr^ 



r-dv 

V h— + cw 
dr 



0, 



(23) 



We adopt the equation of state for the ideal gas: P = poTZT, Ei — 3/2 71T, where TZ — k/nip 

p + pjc^ 

is the gas constant, L is the luminosity, h{T) — — ~ 1 since for the conditions typical 

Po 

for any line-driven wind, quantities P/c^, Eipa/c? are vanishingly small and p — po+0{v^/c'^). 
Furthermore we assume that the gas is isothermal and hence: 



P = aVo, 



(24) 
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2 fdP' 
where a = — — 

\dpo/ T 

(24) to transform the dP/dr term, after some manipulation, from (23) we obtain: 



TZT. Using the continuity equation (10) and the equation of state 



v: \ dv 
1- ^]v— 

dr 



:k 



y/h^nr^ \aeVthM 



GM / 1 
r'^\/h \\/h 

471 Y 



r 1 - ( - 

r 



r-dv 

V ll- h cw 

dr 



w \ 

vn) 

a 



0, 



(25) 



where 



dP 



TZT 



~ and r = Lae/A'KcGM. 



Analogously to CAK theory, equation (25) has a critical point that is not a sonic point, 
which is evident from the fact that, even li v — Vg the last term containing dv/dr does not 
vanish. Thus we are looking for a solution that starts subsonically from r^, goes smoothly 
through the critical point and then reaches terminal velocity v°° at infinity. Converting 
equation (25) to nondimmensional units according to the following formulas: 



X 



V 

Vr 



(26) 



and for simplicity omitting tilde, equation (25) reads: 



f(p.«.-)-(i-^)«J' + ^(;^-r 



X- 



'Vh 



'Vh 



P 



Vh + 



2Vhx^ 



0. 



2 1 
X 2hx'^ 



(27) 



where p — dv/dx, C = c/vc, Oi = a/vc, and /i 



1^/8^ 
2 \Mavth 



which determines the rate 



of an outflow. Equation ( 27) is nonlinear with respect to p, so that a special technique must 
be used. This treatment is well known and had been used in CAK theory. Equation (27) 
may have zero, one or two roots depending on various arguments that are input into it. The 
position of the critical point is determined by the "singularity condition": 



(28) 



The velocity gradient must be continuous in the whole domain of interest thus requiring, 
that the second derivative is defined in the critical point. The regularity condition reads: 
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( 




(29) 



For a given position of the critical point = r^/rg there are three parameters which should 
be determined: fi, p^, v^- The equation for the third parameter is obtained from (27) 
(calculated at the critical point Xc): 



Obtaining of these equations is straightforward but tedious and we derive them in Appendix. 



The procedure of solving of equations (27)- (30) is straightforward. For a given position 
of the critical point Tc we can calculate the value of the velocity and velocity gradient in 
the critical point. Equations (28)- (30) are used to calculate Pc and = Vg/ai. Since 
C = {c/vs)ai, there is only one independent parameter Oi. Then adjusting the position of 
the critical point we integrate the equation (27) inward, looking for the solution that satisfies 
the inner boundary condition. 

To compare self-consistently SLDW and GEF solutions, they should both be matched 
to a flow at rest at some given radius rj„. When a solution for a stellar wind is calculated, a 
photospheric boundary condition is usually adopted. In such a case the position of the critical 
point is adjusted in order to obtain a solution that gives the position of the photosphere 
^ph = ^{'T — 2/3) at a given radius R, which is identified with the radius of a star (see, 
e.g. Bisnovatyi-Kogan 2001). However, in case of disk-powered winds this procedure is 
clearly non-physical and we should approach a different strategy. It is not possible to fit self- 
consistently a solution for the spherically-symmetric wind with that of the accretion disk. 
Generally, it would require a 2D modelling which is beyond of the current studies. 

To be able to compare self-consistently the GEF solution with SLDW solution we should 
start integration in deeply subsonic {v <^ Vg) region, from some initial density This is 
equivalent to the problem of the continuous fitting of the wind solution with that of a static 
core when calculating a structure of a star when mass loss is taken into account. When a 
solution for a stationary outflowing stellar envelope is fitted to that of a static core only T 
and p should be matched at a fitting point rj„: T'^^'"{rin) = T'^°^'^{rin), p'^^'"{rin) = p'^°^'^{rin) 
(see, e.g. Bisnovatyi-Kogan & Dorodnitsyn 1999), where the the fitting point is located in a 
deep subsonic region v{rin) -C Vg- 



0. 



(30) 



3. 



Wind structure 
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We calculate the wind solution in two stages. First we find the subcritical part of the 
solution. On that stage the two-boundary value problem must be solved in order to give the 
position of the critical point, that is adjusted in such way that the solution satisfies the inner 
boundary condition. For the inner boundary condition we prescribe pj„ at rj„, provided, that 
v{rin) <^ Vs- After some experimentation we found that the relative error is indeed becoming 
vanishingly small as the velocity reduces to Vin < 0.01 w^. For the subcritical part of the 
integration domain we have: [xin,x = x^ — 8x], [v{rin)iV — 1 — \pc\^x], and Xin = rin/vg, 
provided that Sx is small enough. For a given position of the critical point, Pc is For a given 
position of the critical point Tc, Pc and Vc are calculated from (36) and (37) of Appendix. 

The critical tool of our calculations is the relaxation method with automatic allocation 
of mesh points. We found that a standard approach based on Runge-Kutta solvers is not 
appropriate in our case because it is very difficult to obtain the desired accuracy when fitting 
to the inner boundary condition using the shooting strategy. There exist a wide variety of 
relaxation methods for the solution of BVP, for example the method of Heney is widely used 
in stellar evolution calculation. Our original code is based on the prescriptions of Press et al. 
(1992) and allows to obtain the solution that satisfies the boundary condition simultaneously 
adjusting the position of the critical point. The position of the critical point x^. is treated 
as an additional variable, as described in Press et al. (1992), 300 grid points are used in the 
subcritical domain. 

After the subcritical part of the solution is fixed we step from the critical point outward 
{v = 1 + \pc\ Sx, X = Xc + Sx), and integrate equation (27) to large radii obtaining v°°. 

The results of the numerical integrations are present on Fig.l and in the Table 1, where 
the following notations are adopted: A^^^ = a;^"^^ - Xc,gef, A^^^^ = xf^^"^ - Xc,gef; 

^MLDW = (^GEF - ^MLDw)/^MLDW ^nd Ag^^ = (^GEF " '^CAk)/^CAK- Each SCt of SolutioUS 

(sj) is characterized by the position of the GEF critical point. Results on Fig.l were obtained 
for a case of a = 1/2. This case is especially suitable for integration because in that case 
the equation of motion can be transformed to quadratic equation with respect to v'. Curves 
are calculated for the following set of parameters: F = 0.5, k — 0.03, Mbh — IO^Mq, 
T = 4 • lO^X. 



3.1. Standeird Line-Driven Wind solution 

The properties of the CAK solution is studied in great detail and it is illustrative to 
present some of these results in the notations adopted in this paper. Thus, the equation of 
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motion reads: 



F(p,.,x)=(l-^).p+^(l-r)-a?^ 



-S^^^^ =0- (31) 



In the SLDW case the velocity gradient in the critical point can be expressed explicitly 
^— ^ = l/xc- The value of Oi = Vg/vc is found from the following relation: 

where e = fs/c, and all other quantaties have the same meaning as in (27). After pc, and 
ai have been found the value of is calculated from the following: 

,,CAK _ (1 — af) {Pc^^Y " /oq\ 
^ - «^2(l+a) ■ (3dj 

As it clear from the Fig.l, the SLDW critical point is situated considerably farther down- 
stream then the GEF critical point, moreover in that part of the domain, where GEF is 
important, v'^^^ « v'^^^. For example, for x^^^ = 50, x'^^^ = 75 and v^^^/v'^^^ ~ 37 
(c.f., Dl). The velocity in the critical point spans the following range of values: vjvth — 1541 
for Xc — 22.5 drops to Vc/vth — 674 for Xc ^ 150. 



3.2. Modified Line-Driven Wind solution 

It is important to understand what is the relative impact to the total dynamics, of the 
acceleration that is due Sobolev effect ( when only Doppler effect is taken into account) 
alone in comparison with GEF case. In other words we would like to see what happens 
if we take into account only Doppler effect, neglecting gravitational redshift, but treating 
the rest of the problem in General Relativity. As it was emphasized in the Introduction 
this approach is not self-consistent, nevertheless considering such a solution can give us an 
important insight to the relative importance of the effect. In such a case the equation of 
motion will be identical to equation (27) apart of the last term, which in such a case reads: 

,MLDW f VX^ 



In Dl it was found that using the Paczynski-Wiita (PW) modified potential results in an 
increase of the terminal wind velocity. It was considered to be natural, because a wind needs 
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to have steeper gradient in order to move out of the sharper potential well. As it is clear 
from Fig. 1, making use of general relativity increases this effect (c.f. Dl, Fig. 4). However, 
for a given radius, the position of Modified Line-Driven Wind (hereafter MLDW) critical 
point differs only slightly from that of GEF: A^^^dw ^ 2,mldw _ ^gef _ q q2 ^ o.7 for a 
considered range of x^'^^ . Note, that for Xc = 15 ^ 100, critical velocity of MLDW is almost 
four times greater then in the GEF case. For example, for Xc — 15, Vc/vth — 99, and Xc — 100 
Vc/vth — 98. Thus our recent studies confirms the qualitative picture which was obtained in 
Dl. 



3.3. Gravitationally Exposed Flow 

In the general rclativistic calculation presented here wc find a considerable gain in the 
wind terminal velocity in comparison with both the CAK solution and the semi-classical 
solution of Dl. 



Model 




A CAK 
c 


A MLDW 
c 


A OO 

^MLDW 


A OO 

^CAK 


Sl 


15 


7.47 


0.016 


0.23 


0.57 


S2 


20 


9.95 


0.03 


0.232 


0.51 


S3 


50 


25 


0.16 


0.24 


0.37 


S4 


100 


50 


0.66 


0.29 


0.2 



Table 1: Comparison of GEF solution and SLDW solution. See text for details. 

The position of the GEF critical point is found to be closer to BH than the SLDW 
critical point. From the Table 1. we see that even when the solution originates (in fact 
it is determined by the position of the critical point) sufficiently far from BH there exist a 
valuable gain in v*^^^ in comparison with CAK case. 

4. Discussion and conclusions 

Line-driven winds represent the most plausible explanations of fast outflows from various 
astrophysical sources. 

The most successful theory presented so far has been that of Castor, Abbott and Klein 
(1975), (CAK), where a theory of an 0-type star wind was developed. A modification of this 
approach was used to explain outfiows in active galactic nuclei. It is now widely accepted 
that most of AGN manifests itself by fast uncoUimated outfiows. Such fast (up to 0.1c) 
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winds are believed to be driven by the radiation pressure on spectral lines. An accretion 
disk around supermassive black hole (BH) is believed to be a source of such a powerful 
continuum radiation. A lot of work have been done to put together this jig-saw puzzle (see, 
e.g. Arav & Li (1994); Arav, Li & Begelman, (1994); Murray et al. (1995); Proga, Stone & 
Drew (1998); Proga, Stone & Kallman (2000) and Stevens & Kallman (1990), note that this 
problem is closely connected with that of modelling of winds form massive X-ray binaries) 

It is in the paper of Dl) where the mechanism of plasma acceleration due to absorbtion 
of the radiation flux in lines in a strong gravitaional field was first investigated. A photon 
emitted deeply in the potential well will suffer a continuous redshifting of frequency, that 
should be additionally taken into account if a wind is accelerated near compact object. A 
parcel of gas sees matter that is upstream as being redshifted both due to difference in 
velocities (as in classical Sobolev approximation) and due to gravitational redshifting which 
exist regardless on whether gas is moving or not. As it had been demonstrated in Dl, the 
radiation pressure on spectral lines becomes a function of the local velocity gradient and 
the gradient of the gravitational potential. Thus it has been concluded that the greater 
the gravitational redshifting, the more effectively a line is shifted to the extent where the 
radiation flux is un-attenuated by the line opacity. Since in such a case the gravitational field 
'exposes' a wind to the un-attenuated continuum, we call this kind of flow 'gravitationally 
exposed flow' (GEF). 

The generalization of these studies in the frame of General Relativity (GR) is the prob- 
lem that has been addressed in this paper. Only this approach allows to take self-consistently 
into account both Doppler and gravitational redshifting. The main goal of these studies was 
to confront GEF regime with GAK wind. 

Using the Sobolev approximation, generalized within the framework of GR, the acceleration 
that is due to absorbtion in a single line was found to be 



9i ~ iVdoQ-; — I" c— — ) that should be compared with the CAK case: ~ — . The 



acceleration due to gravitational redshifting is most important at the bottom of the wind 
where the velocity is small. In our treatment terms of the order 0(u^/c^) had been neglected, 
thus it is roughly accurate for velocities as fast as ~ 0.2c. However the relative lost of 
accuracy at mildly relativistic velocities is not very important, because our goal was to 
calculate the relative importance of the effect (which is important at low v). The derived 
general relativistic equation of motion has a critical point that is different from that of 
CAK (Note that the CAK point is not a sonic point). Since this equation is nonlinear with 
respect to velocity gradient the CAK approach to such an equation was adopted. In order 
to compare GEF solution with the SLDW we numerically solve the two boundary value 
problem. A relaxation method is found to be a very important tool to flnd such a solution. 
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An integral impact of the gravitational redshifting to the radiation acceleration can give 
a considerable gain in terminal velocity (Av°°/v°° > 0.4). We defer a detailed study of 
the mathematical properties (including the analysis of stability) of the GEF equations to a 
separate paper. 

Finally we summarize most important results which have been obtained in the current 
studies: 

1. A wind driven by the radiation pressure on spectral lines was considered in the frame of 
General Relativity. Following Dorodnitsyn (2003), we argue that it is important to take 
into account the gravitational redshifting of the photon's frequency, when calculating 
the radiation force. 

2. A generalization of the Sobolev approximation in GR was developed and the general 
relativistic equation of motion with the radiation pressure force on spectral lines was 



3. The results of the numerical integration of the equation of motion demonstrate that 
taking into account gravitational redshifting can result in a wind that is considerably 
more fast than previously assumed on the ground of the GAK theory. 
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Appendix 



From the 'singularity' condition (28) we obtain the value of /i: 




(35) 



where he = h{xc) Substituting (35) to (30) and solving the resultant equation with respect 
to Pc will result: 



_ -ai^ {a + AhcXca- () +( (-1 + ct C - C T) 

2 {ai^ -1) hcXc^ (a-l) ■ ^ ^ 

After a lengthy calculus of the 'regularity condition' (29), taking into account (35) the 
quadratic equation can be obtained: 



C2Pc + ClPc + Co = 0, (37) 
where q should be calculated from the following: 



Co = - [ai^ [Shc^ a -2hica + hic C)) 

+ hicC ((-2q;C) + v^aCr) 

Cl = 2 {a,^ -1) h, {h,Xc^ {2a-l)+ aC) 
C2 — 8 ai^ b Xc ot 



where hi^. = dh/dx{x^. For a given position of the critical point Tc equations (36) and (37) 

are used to calculate pc = {dv/dx)c and ai — Vs/vc- The mass-loss rate is straightforwardly 

, , , , SttGM fTk^ 
calculated from M = — 
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Fig. 1. — Solutions of the equation of motion (27). Solid line - " Gravitationally Exposed 
Flow" (GEF) solution, dashed line - standard line-driven wind (SLDW) regime. Stars indi- 
cate GEF critical points, circles - SDLW critical points, cf. Figure 1. of Dorodnitsyn (2003) 
Labels Si - S4 mark solutions with different locations of the critical point. The following set 
of parameters was adopted: F = 0.5, k = 0.03, Mbh = lO^M©, T = 4 • WK. 



